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Simulations of collisionless oblique propagating slow shocks have revealed the existence of a tran- 
sition associated with a critical temperature anisotropy e = 1 — /J>o(P\\ — Pi_)/ B 2 — 0.25 (Liu, Drake 
and Swisdak (2011) pQ). An explanation for this phenomenon is proposed here based on anisotropic 
fluid theory, in particular the Anisotropic Derivative Nonlinear-Schrodinger-Burgers equation, with 
an intuitive model of the energy closure for the downstream counter-streaming ions. The anisotropy 
value of 0.25 is significant because it is closely related to the degeneracy point of the slow and inter- 
mediate modes, and corresponds to the lower bound of the coplanar to non-coplanar transition that 
occurs inside a compound slow shock (SS) /rotational discontinuity (RD) wave. This work implies 
that it is a pair of compound SS/RD waves that bound the outflows in magnetic reconnection, in- 
stead of a pair of switch-off slow shocks as in Petschek's model. This fact might explain the rareness 
of in-situ observations of Petschek-reconnection-associated switch-off slow shocks. 

I. INTRODUCTION 



Shocks in isotropic MHD have been intensively studied, including the existence of intermediate shocks (IS) (EHS 
and references therein) , the occurrence of dispersive wavetrains [BHB] , and the nested subshocks inside shocks predicted 
by the Rankine-Hugoniot jump conditions [9[|T0]. In a collisionless plasma, the effects of temperature anisotropy need 
to be considered, which can be done for linear waves with the Chew-Goldberger-Low (CGL) framework [TTJ [12]. Hau 
and Sonnerup have pointed out the abnormal properties of the linear slow mode under the influence of a firehose-sense 
(P|| > P_\_) pressure anisotropy, including a faster phase speed compared to the intermediate mode, a fast-mode-like 
positive correlation between magnetic field and density, and the steepening of the slow expansion wave. In kinetic 
theory both anisotropy and high /? can greatly alter the linear mode behavior p~3j [14] . The anisotropic Rankine- 
Hugoniot jump conditions have been explored while taking the downstream anisotropy as a free parameter [T5HT8] . 
while Hudson [17] calculated the possible anisotropy jumps across an anisotropic rotational discontinuity. Karimabadi 
et al. [14] noticed the existence of a slow shock whose upstream and downstream are both super-intermediate. But, 
a comprehensive nonlinear theory describing the coupling between slow and intermediate shocks under the influence 
of temperature anisotropy has not yet been presented. 

In Petschek's description of magnetic reconnection, the reconnection exhaust is bounded by a pair of back-to-back 
standing switch-off slow shocks. Particle-in-cell (PIC) simulations of such shocks [TJ [19] exhibit large downstream 
temperature anisotropics. In Liu et al. (2011) pQ (hereafter called Paper I) and Fig. [I] of this paper we show that 
when the parameter e = 1 — /J>o(P\\ — P±)/B 2 = 0.25, the behavior of the coplanar shock undergoes a transition to 
non-coplanar rotation. This firehose-sense temperature anisotropy slows the linear intermediate mode and speeds up 
the linear slow mode enough so that, at some point, their relative velocities can be reversed [TTJ [12]. This reversal 
is reflected in the structure of the Sagdeev potential (also called the pseudo-potential) [20], which characterizes the 
nonlinearity of the system. In this work a simplified theoretical model is developed to explore the effect of temperature 
anisotropy on the structure of the Sagdeev potential and to provide an explanation for the extra transition inside the 
switch-off slow shock (SSS) predicted by isotropic MHD. The theory suggests that in PIC simulations a compound 
slow shock (SS)/ rotational discontinuity (RD) is formed instead of a switch-off slow shock. This work may help to 
explain satellite observations of compound SS/RD waves [21, 22 , anomalous slow shocks [23 j and the trapping of an 
RD by the internal temperature anisotropy of a slow shock in hybrid simulations [24] . 

In Sec. II of this paper we introduce our model equations for studying the nonlinear coupling of slow and interme- 
diate waves under the influence of a temperature anisotropy. In Sec. Ill we calculate the speeds and the eigenmodes 
of slow and intermediate waves. Sec. IV points out the existence of extra degeneracy points between slow and in- 
termediate modes introduced by the temperature anisotropy, and comments on the consequences (in the context of 
the Riemann problem) of having the slow wave faster than the intermediate wave. In Sec. V we introduce a simple 
energy closure. In Sec. VI. A. we calculate the pseudo-potential of stationary solutions, and apply the equal-area rule 
to identify the existence of compound SS/RD waves and compound SS/IS waves. In Sec. VI. B we demonstrate the 
significance of e = 0.25 as being the lower bound of the SS to RD transition in compound SS/RD waves. In Sec. VII 
we discuss the time-dependent dynamics that help keep e = 0.25. In Sec VIII we provide more evidence from PIC 
simulations to support the existence of compound SS/RD waves at the boundaries of reconnection exhausts. In Sec 
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IX, we summarize the results and point out the relation between compound SS/RD waves and anisotropic rotational 
discontinuities [T7] . 

II. THE ANISOTROPIC DERIVATIVE NONLINEAR SCHRODINGER-BURGERS EQUATION 

Instead of analyzing the anisotropic MHD equations, which have seven characteristics (waves), we simplify the 
system into a model equation that possesses only two characteristics. This model equation will be ideal for demon- 
strating the underlying coupling between the nonlinear slow and intermediate modes. Beginning with the anisotropic 
MHD equations [15] , we follow the procedure of Kennel et al. [4j [25 j to derive the Anisotropic Derivative Nonlinear 
Schrodinger-Burgers equation (ADNLSB) (see Appendix I for details), 
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This equation describes waves that propagate in the ^-direction in the upstream intermediate speed frame. In this 
frame 77 = x' — ^/EoCAnt is the spatial coordinate with C\ n = ^/(/ioPo), and r = ^/EoCAnt is the time used to measure 
the slow variations (such as steepening processes), = ~B t /B x , where the subscript "f represents the component 
tangential to the wave- vector and here will be in the y-z plane. The anisotropy parameter e = 1 — fio(P\\ — P±)/B 2 . 
The subscript "0" denotes the upstream parameters. The right hand side term proportional to the ion inertial length, 
di = <s/ 'mi j (/lone 2 ) , represents dispersion (which can be viewed as the spreading tendency of Fourier decomposed 
waves of different wavenumbers), while the term containing R describes dissipation from magnetic resistivity. Here R 
is a constant. The terms proportional to a and ft are the nonlinearities of this wave equation, where 
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and 
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Here C| = 7P0/P0 an d 7 = 5/3 for monatomic plasma. Since we are studying reconnection exhausts, #0 (the angle 
between the upstream magnetic field and e^) is typically large (~ 80°). Therefore f3 n = f3/ cos 2 0q = C\jC 2 An ^> 1, 
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-Cg < 0. Hence this equation describes only the slow and intermediate modes [4 (this 



is shown explicitly in the next section). This fact relates to the degeneracy properties of ideal MHD for parallel 
propagating waves, namely that the fast and intermediate modes degenerate in j3 < 1 plasmas, while the slow and 
intermediate mode degenerate in j3 > 1 plasmas. Finally, we note that Eq. ([T]) is applicable in the weak nonlinearity 
limit. 



III. THE CONSERVATIVE FORM- WAVE PROPAGATION 



In order to explore the structure of the reconnection exhaust, a comprehensive understanding of how waves connect 
to each other across a transition is required. This is called a Riemann problem. Neglecting the source terms on the 
right hand side (RHS), the left hand side (LHS) of Eq. ([!]) is a hyperbolic equation in conservative form. 

Letting e — So = Se(b z , b y ), and b^o = b z oe z , we obtain, 
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We can obtain the characteristics (waves) of this equation by analyzing its flux function, f. Its Jacobian is 
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where 8e bz = d(5s)/db z , 5e by = d{5e)/db y . One eigenvalue (also called the characteristic speed) is 

Asl = a(3b 2 t - b 2 t0 ) + Q(Se + b z 5e bz + b y 5e by ) 

with eigenvector 
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In isotropic ideal MHD, where the (e — Sq) term is dropped, the eigenvalue in the infinitesimal limit (b t —> b t o) is 
Xsl = 2a6f , which is the phase speed of the linear slow mode in the intermediate mode frame. The subscript "SL" 
means the slow mode. The eigenvector indicates that slow mode is coplanar (i.e., in the radial direction in b z — b y 
space). 

The other mode has eigenvalue 



A/ = a(b 2 t - b 2 t0 ) + nSe 



and eigenvector 
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In isotropic ideal MHD, where the (e — eo) term is dropped, the eigenvalue in the infinitesimal limit is A/ = 0, 
which is the phase speed of the linear intermediate mode in the intermediate mode frame. The subscript "I" means 
the intermediate mode. The eigenvector indicates that this intermediate mode is non-coplanar (i.e., in a non-radial 
direction) . 

It can be shown that V q A/(q) • 17 (q) = for all q. This means that along the eigen-direction of the intermediate 
mode the characteristic speed is constant, and thus the mode exhibits no steepening or spreading, just as is the case for 
its counterpart in isotropic MHD. (This behavior is also confirmed by the anisotropic MHD simple wave calculation; 
see Appendix II) . Therefore the intermediate mode is termed "linearly degenerate" . If we are looking for a transition 
in the — r^x direction (toward b t = 0), then the portion of the slow mode with V q Asx(q) • i*sx(q) < will steepen 
into a slow shock. When Xsl at the downstream of a transition is larger than that at the upstream, the downstream 
wave will catch up with the upstream wave and thus steepen. 



IV. A NEW DEGENERACY POINT DUE TO THE TEMPERATURE ANISOTROPY 



In the Riemann problem for our two mode system, we seek to determine the middle state q m that connects the 
faster "2- wave" from a given state q r , to the slower "l-wave" from a given state qi (see Fig. [2] (a); the subscripts "r" 
and "1" mean right and left respectively). In order to determine the path that connects q r to qj in the state space 
(b z — by space, in this case), the Hugoniot locus that connects q^ or q r to a possible asymptotic state by shock waves 
needs to be calculated, as do the integral curves for possible rarefaction waves (see, for example, [26]). The Hugoniot 
locus in state space is a curve formed by allowing one of the parameters in the standard Rankine-Hugoniot jump 
condition to vary. The integral curve is formed by following the eigenvector from a given state in state space. 

In order to proceed we further assume a gyrotropic energy closure, 5e = Se(b t ), which allows us to write 
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The Hugoniot locus and integral curves of the intermediate mode are identical in our system. This is also the case 
for the slow mode (we perform the calculation in Appendix III). For a given state, q r = h t Q = (fr^O), the Hugoniot 
locus and integral curve of the intermediate mode are 
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which is a circle in state space. Note that even though we can calculate the Hugoniot locus and integral curve for the 
intermediate mode, the solution is the same as for a finite amplitude intermediate mode that does not steepen into a 
shock or spread into a rarefaction. For the slow mode, the Hugoniot locus and integral curve are 

by = (14) 

which is in the radial direction in state space. This direction implies that the slow shock is coplanar, even in the 
presence of temperature anisotropy, just as is the case for its counterpart in full anisotropic MHD [15]. In the isotropic 
case, this curve forms a slow shock if the path is toward the origin, and a slow rarefaction if the path is away from 
the origin. 

The state q r can connect to q/ = (0,0) by following the Hugoniot locus of a slow mode that starts from q r , as 
shown in Fig. [2^b). In isotropic MHD this forms a switch-off slow shock. However, a strong enough temperature 
anisotropy introduces new degeneracy points (which occur where Xsl — A/ = 0) when 2ab t + QSe^ = 0, other than 
the traditional degeneracy point at b t = 0. These points form a band circling the origin as shown in Fig. [2jc). Inside 
the band, the intermediate mode is slower than the slow mode. Physically, this implies that a rotational intermediate 
mode can arise downstream of a slow mode, something which is not allowed in a Riemann problem in isotropic MHD. 
This effect is realized when the path along the Hugoniot locus (— tsl direction) of the slow mode from q r switches to 
the solution of the intermediate mode (circular direction) somewhere (q m ) inside the degeneracy band. 

This behavior can explain the morphological differences between the shock simulations in the cases Obn = 30°, 45° 
and those for Obn = 60°, 75°, 83° of Paper I [1 . The latter has an extra transition to the rotational direction that 
is similar to the path in Fig. [2|c). We now look for a similar effect in state space and a way of determining q m in a 
more detailed analytical model. 



V. AN ENERGY CLOSURE BASED ON COUNTER-STREAMING IONS 



In order to close the ADNLSB equations, we need a energy closure e(b t ). The modeling of the energy closure for a 
collisionless plasma has historically been difficult. The Chew-Goldberger-Low (CGL) condition [27] is one choice, but 
it does not work well when streaming ions are present. Since we are here just trying to qualitatively demonstrate the 
underlying physics, we will assume that we have a e(b t ), where e and \B\ are simply related by 

e = ci - (15) 

with positive constants c\ and C2 and the condition c\ — C2/B 2 = so is imposed. This functional form is motivated 



by the nearly constant parallel pressure maintained by free-streaming ions (02 ~ P\\). Although Eq. (15) is strictly 
empirical, results from PIC simulations (see Fig. [3| suggest that C2 = 0.5 provides a reasonable first approximation 
and will be used in the following calculations. 
We take the variation, 

fe = £ ~ £o ~ ^ m = W^ l ~ hl) (16) 

where b 2 = 1 + b 2 . This parameterization will be valid whenever Se(b t ) <C £o- Therefore, an effective nonlinearity in 
Eq. ([!]) can be written as, 

a ef ?(b t ) = a + n-^[ (17) 

The most important conclusions in the remainder of this work do not depend on the details of the closure, but only 
that e decreases as \B\ decreases. Note that Q ~ l/2eo is mostly positive in the limit in which we are interested. This 
fact will be used in the following section. 



VI. THE PSEUDO-POTENTIAL: LOOKING FOR A STATIONARY SOLUTION 
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A. The formation of compound SS/RD waves and SS/IS waves 

In order to determine both where the path in the state space of Fig. [2|c) will turn to the intermediate rotation and 
the nontrivial coupling of the slow and intermediate modes when temperature anisotropics are present, we construct 
the pseudo-potential of a stationary solution. We will look for a equation that possesses traveling stationary waves, by 
substituting = h t [£(rj — Vst)] (where Vs is the speed of the stationary wave observed in the upstream intermediate 
frame) into Eq. ([!]) and integrating over £ once. We obtain 

Rdzbt - -^—did^ex x b t ) = -V s (b t - b t0 ) + a^(b t )b t (b 2 - b 2 t0 ) = -F = d ht V. (18) 



In this formulation we can treat bt as a spatial coordinate, and £ as time. The terms on the LHS of Eq. (18) behave 
analogously to, respectively, a frictional force and a Coriolis force with rotational frequency di/^/Eo and rotational 
axis e x . A pseudo-potential \I/ that characterizes the nonlinearity is uniquely defined because dF z /db y — dF y /db z . 
We are only interested in the small Vs limit, because the upstream values with subscript "0" are expected to be the 
upstream values of a switch-off slow shock in ideal isotropic MHD, which propagates at the upstream intermediate 
speed. The anisotropy in our PIC simulations does not seem to significantly change this behavior pQ. 



Calculating the pseudo-work done on the pseudo-particle, J[Eq.(18)] • d^h t d^ we obtain 



ydown 

*lup wn = R / (W 2 ^ < (19) 

J up 

Note that from upstream to downstream is in the negative £ direction. The pseudo-particle will move to a lower 
potential, while its total energy is dissipated by the resistivity and the rate of the drop depends on the strength of the 
resistivity. Kennel et al. [4! have shown that when pseudo-particles move toward lower pseudo-potentials, the entropy 
increases and so the resulting shock is admissible. Note that the Coriolis-like force does not do work. It only drives 
rotation of the pseudo-particle on the iso-surface of the pseudo-potential and hence causes stable nodes to become 
stable spiral nodes and unstable nodes to become unstable spiral nodes, thus leading to the formation of dispersive 
wavetrains [TJ [8] . We will neglect its effect in the following discussion. 

The pseudo-potential is shown in Fig. Qa) for the parameters 6$ = 42°, /3 = 1, £o = 1> c 2 = 0.5 and Vs = 0. The 
temperature anisotropy has turned the origin from a local minimum of the pseudo-potential in the isotropic MHD 
model to a local maximum. We term this the "reversal behavior". A pseudo-particle initially at point q r will slide 
down the hill in the slow mode eigen-direction, and then follow the circular valley in the intermediate eigen-direction. 
Without the reversal behavior (e.g., in isotropic MHD) the pseudo-particle will slide down to the origin and form a 



switch-off slow shock. The trajectory of the pseudo-particle can also be calculated by numerically integrating Eq. (18) 
with respect to £. In (b) the variation of the temperature anisotropy is shown. Similar reversal behaviors can be 
found in fully anisotropic MHD with the energy closure used here, or with the CGL closure (see Appendix IV. A). 

In Fig. [4] (c), we plot a cut of the pseudo-force, F z , effective a e R and the pseudo-potential \I/ along the b z axis 
(by = 0), which is the eigen-direction of a slow mode beginning at q r = (b z o,0). Here 

F z = -a eE (b z )b z (b 2 z - b 2 z0 ) + V s (b z - b x0 ). (20) 

It is clear that ^ m in occurs at F z = in Fig. [4jc), since \I/ is constructed by integrating the pseudo-force F. In Fig. [4] 
(d) we plot cuts of the characteristics of the slow and intermediate modes along the b z axis. 



XsL = a(3b 2 z -b 2 z0 ) + Q^ 



(21) 



A/ = a eS (b z )(b 2 z - b 2 z0 ) (22) 

The temperature anisotropy has changed the structure of these characteristics. As a result, there are new degeneracy 
points (Xsl = Xi) between slow and intermediate waves such as the point "D". The slow characteristic shows extra 
nonconvexity points, where no steepening and spreading occurs (i.e., V q Asx(q) • i*sx(q) = 0), such as the point a 
(the local maximum of Xsl)- This is clearer when we compare the slow characteristic here to that in the isotropic 
case shown in Fig. [lOj where b t = is the only degeneracy point and the nonconvexity point of the slow mode. 

In order to identify the nonlinear waves determined by the route of the pseudo-particle, we apply the equal- area 
rule, which tells how shocks are steepened from characteristics. The equal-area rule (see Appendix V for more details) 
applied to Xsl shows that the sliding route (point q r to q m ) forms a slow shock. Since V q Asx(q) • i*sx(q) < 
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therefore Asx|q m > Asx| qr , and thus the slow mode will steepen until the red area above the horizontal line Vs = 
equals the red area below Vs = 0. The slow shock transition immediately connects to the intermediate mode (point 
q m to b, which is also from the equal-area rule on Xsl) in the valley. The fact that both the upstream (point q m ) 
and downstream (point b) travel at the local A/ makes the intermediate discontinuity a RD. By comparing (c) and 
(d), we note that the potential minimum is exactly the location of q m as expected and it is below the degeneracy 
point (b z \o > &z|q m )- This fact is consistent with the comment in section V, which predicts that q m will be inside 
the degeneracy band. The horizontal lines q r — q m and q m — b measure the propagation speed of the SS and the 
RD, which in this case are both zero in the upstream intermediate frame. They therefore form a compound SS/RD 
wave. The downstream of the slow shock (point q m ) is not able to connect to the slow rarefaction (SR) wave (point 
a) and thus not able to form a compound SS/SR, since the rarefaction is faster than the shock itself. This model gives 
an theoretical explanation for the possible satellite observations of compound SS/RD J2TJ [22] , and the "compound 
SS/RD/SS waves" seen in hybrid simulations [24] . 

When Vs > the potential tilts down in the negative b z direction (see Fig. [5]). In this case, the slow shock (point 
q r to q m in Fig. [5^d)) with shock speed Vs is connected by an intermediate shock (IS) (point q m to b) with shock 
speed V5, whose upstream is super-intermediate (Vs > Aj| qm ) while the downstream is sub-intermediate (Vs < A/|b). 
This forms a compound SS/IS wave. Note that the intermediate shock is not steepened from the intermediate mode 
(which is consistent with the discussion in Sec. IV), but is steepened from the slow mode. The slow shock is abnormal 
with both upstream (point q r ) and downstream (point q m ) being super-intermediate. Karimabadi et al [14] call this 
kind of slow shock an anomalous slow shock. When Vs < the potential tilts up in the negative b z direction and 
there is no extra transition at the SS downstream, since ^(b) > \£(q m ) in this case and is therefore not accessible. 

These results are independent of the details of e(b t ) : but only require the reversal behavior somewhere downstream 
of the slow shock. This fact can be inferred from a simple relation: Vs — \i(b z ) = (b z o/b z )Vs for b z = ^| qrn or 6 z |b, 
regardless of the detail of s(b t ) (from Eq. (20) and ([22])) . When Vs = 0, this relation ensures that the SS can always 
connect to a RD since Vs — Aj = at both points q m and b. When Vs > 0, the SS can always connect to an IS since 
Vs — A/ is positive (super- intermediate) at q m and negative (sub- intermediate) at b. We therefore conclude that the 
abnormal transitions of magnetic field structures seen in the PIC simulations of Paper I are most likely the transitions 
from the SS to the RD in a compound SS/RD wave or the SS to the IS in an SS/IS wave. We can hardly distinguish 
between these two compound waves in our PIC simulation, since Vs is small and the time-dependent dynamics add 
uncertainties in measuring the exact value. We focus on further analyzing the compound SS/RD wave. 



B. The significance of e — 0.25 

For SS/RD waves (Vs = 0), the stationary points along b ZjS are the roots of F z = 0, 

a eH (b Z:S )b Z:S (bl s -b 2 z0 )=0 (23) 

Here the subscript "s" represents "stationary". We have three traditional stationary points, b ZiS = b z o (point q r ), 
—b z o and 0, as well as a new stationary point due to the temperature anisotropy, b z ^ m (the transition point q m ) 
determined by a e ff(^,m) = 0- The fixed-point analysis of the first three points in isotropic fluid theory can be found 
in literature [10[H][1H]. 

As shown in Fig. [6] there is no slow mode transition if 

3^-4 

so < = 0.25 (24) 

07 — 1 

(with 7 = 5/3 for monatomic plasma; note that this relation is independent of Obn and /3), which occurs when the 
nonlinearity a of Eq. changes sign from a < to a > 0. A positive a will result in a positive a e ^ in Eq. (17), 



and therefore no solution for b Z:Tn . Only rotation of the magnetic field is thus allowed. If Eq > 0.25, we can further 
show that Eq > e m (= £|q m ) > 0.25 is always true for a slow shock transition from eq to e m in this compound wave 
by the full jump conditions of anisotropic-MHD (Appendix IV. B). Therefore, the nonlinear fluid theory provides a 
lower bound of £ m > 0.25 at the SS to RD transition inside these compound waves, regardless of the details of e(b t ). 
In other words, the downstream magnetic field cannot exhibit switch-off behavior if the firehose-sense temperature 
anisotropy is strong enough. This fact explains the non-switch-off slow shocks often seen in kinetic simulations [29] 
and satellite crossings [30 . Once it transitions to the intermediate mode, a gyrotropic s(b t ) will stay close to e m , 
since the intermediate-rotation nearly preserves the magnitude of the B field and therefore e. Note the assumption of 
gyrotropic e(b t ) is expected to be valid only in length scale larger than local ion inertial length and ion gyro-radius. 

In these demonstrations we use shocks with moderate parameters, such as Obn = 42° and /3q = 1. In general, larger 
bn 1 Po, and smaller eo will make the ratio |fi|/|a| larger, and therefore generate a stronger reversal tendency. An 
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analysis with full anisotropic MHD (Appendix IV) should be used for strong slow shock transitions, due to the limits 
of the ADNLSB, although the underlying physical picture will be similar. 



VII. TOWARD THE CRITICAL e = 0.25: TIME-DEPENDENT DYNAMICS 



The initial conditions that characterize the exhaust of anti-parallel reconnection (initial b y = 0) require h t = at 
the symmetry line at later time. This eventually forces the pseudo-particle to climb up the potential hill to the local 
maximum (b z = 0,b y = 0), which implies an intrinsic time-dependent process at the symmetry line since Eq. (18) 
does not yield such a solution. Meanwhile, the fact that needs to go to zero at the symmetry line provides a 
spatial modulation on the amplitude of the rotational intermediate mode. Note that the transition point from SS 
to RD in compound SS/RD waves could potentially induce modulation too. As suggested in Paper I, a spatially 
modulated rotational wave tends to break into c^-scale dispersive waves, which can make the rotational component 
of the transition very turbulent. 

As pointed out in Sec VI. B, the nonlinear fluid theory of the time-independent stationary solutions only provides 
a lower bound £ m > 0.25 for the transition point inside these compound SS/RD waves. Counter-streaming ions, 
by raising Py, push e m lower. Once £ m is lower than 0.25, the magnetic field rotates, generates descale waves, and 
scatters Pu into P±. This raises £ m , changing the functional form of 5s and driving it toward 0, which self-consistently 
results in a transition at the potential minimum where a e ff = a = 0, and thus e = 0.25. This argument explains the 
e — 0.25 plateau observed in the PIC simulations for different shock angles (see Fig. [I]). With Se = 0, this point is 
exactly the degenerate point of the slow and intermediate modes. 



VIII. THE SUPPORTING EVIDENCE FROM NUMERICAL EXPERIMENTS 



The evidence for a slow mode connecting to rotational waves can be seen in the PIC simulations that are discussed 
in detail in Paper I. In previous kinetic simulations, the downstream rotational waves were often identified as slow 
dispersive wavetrains arising from, for instance, the second term in the RHS of Eq. ([!]). Here we present further 
evidence, in addition to the numerical evidence that e = 0.25, to support the idea that the downstream rotational 
mode is tied to the intermediate mode. Fig. [7] shows the results from three PIC simulations that were designed 
to explore the structure of reconnection exhausts in the normal direction. The format is the same as for Fig. 5 of 
Paper I, with the first column showing e, the second column the magnetic field components, and the third column 
a hodogram of the fields. The dashed curves in the second column (from ideal isotropic MHD [31]) indicate that a 
pair of switch-off slow shocks or a pair of rotational discontinuities will propagate out from the center. All three cases 
show the correlation between e = 0.25 and the transition from coplanar to non-coplanr rotation of the downstream 
magnetic fields. The hodograms are readily comparable to the state space plots such as Fig. [2](c). In Fig. [7^a), the 
downstream region of a slow shock shows a high wavenumber (~ 6^) left-handed (LH) polarized rotational wave, 
which is difficult to distinguish from the predicted downstream ion inertial scale dispersive slow mode wavetrain [6]. 
Fig. (7^b) shows results from a simulation with a larger initial current sheet width and exhibits a longer wavelength 
(~ 30di) LH rotational wave which can be identified as an intermediate mode The intermediate mode breaks into 
smaller ion inertial scale waves, which have been identified as dispersive waves in Paper I. By comparing (a) and 
(b), we note that the downstream primary rotational wave tends to maintain its spatial scale as an intermediate 
mode with non- steepening and non-spreading properties. Another way to distinguish the dispersive behavior from 
the non-dispersive rotation is by including a weak guide field. In Fig. (7^c), the front of the rotational downstream 
wave turns into a well-defined RD when a weak guide field is included. Its amplitude is about the same as that of the 
large amplitude rotational waves in (a). Most importantly, there is a clear slow shock ahead of the RD. Because the 
symmetry of the initial condition is broken by the guide field, the downstream RD does not need to end at b y = 0; 
instead it ends inside the potential valley at h t — (0,b Zjfn ) (see the hodogram of (c)), as expected. 



IX. CONCLUSION AND DISCUSSION 



The existence of compound SS/RD, SS/IS waves arising from firehose-sense and | B | -correlated e (temperature 
anisotropics) are theoretically demonstrated by analyzing the anisotropy-caused reversal of a pseudo-potential. The 
pseudo-potential is known to characterize the nonlinearity of hyperbolic wave equations. Extra degeneracy points 
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between slow and intermediate modes as well as extra non-convexity points in the slow characteristics are introduced 
by the temperature anisotropy. The slow shock portion of a compound SS/IS wave is an anomalous slow shock with 
both up and downstream being super- intermediate. The nonlinear fluid theory provides a lower bound of e = 0.25 for 
the SS to RD transition, regardless of the details of the energy closure e(b t ). The wave generation from the rotational 
intermediate mode discussed here and in Paper I helps keep e = 0.25. This explains the critical anisotropy plateau 
observed in the oblique slow shock PIC simulations documented in Paper I. This study also suggests that it is a 
pair of compound SS/RD waves that bound the antiparallel reconnection outflow, instead of a pair of switch-off slow 
shocks as in Petschek's reconnection model. This fact explains the in-situ observations of non- switch- off slow shocks 
in magnetotail [30 . It also provides a theoretical explanation of the observations of "Double Discontinuity" by Whang 
et at. [2TJ [22] with GEOTAIL data, and also the step-like slow shocks seen in Lottermoser et aZ.'s large-scale hybrid 
reconnection simulation [29]. In previous hybrid and PIC simulations, the downstream sharp rotational waves were 
often identified as slow dispersive waves of a switch-off slow shock. Instead, we propose that they are the intermediate 
portion of the compound SS/RD wave. The slow shock portion becomes less steep due to the time-of- flight effect of 
backst reaming ions. 

The singularity of e = 0.25 was also noticed by P. D. Hudson [17] in his study on the anisotropic rotational 
discontinuity (A-RD). Unlike the RD in a compound SS/RD wave, an A-RD changes both e and thermal states. 
Through the constraint of the positivity of Py, Pj_ and £? 2 , he derived all of the possible jumps (independent of the 
energy closure used) of the temperature anisotropy across an A-RD, as shown in Fig. |8ja). (Note that in Fig. [8] 
= e up — Sdown — >■ when e — » 0.25.) The ADNLSB inherits most of the hyperbolic properties (such as the 
extra nonconvexity and degeneracy points) of anisotropic MHD and also this singular behavior. We can tell this by 
searching for stationary solutions where the pseudo-force F z = (the general form of Eq. (|20|)). Again, 

V s (b z - b z0 ) - ab z (b 2 z - b 2 z0 ) - Qb z (e - e ) = (25) 

A stationary A-RD exists at Eq = 0.25 {i.e., Vs = 0, a = and b z ^ 0), only if we require e — £o = 0. An arbitrary 
magnetic field magnitude and rotation are then allowed, as shown by Hudson. After further constraining the possible 
jumps by requiring that entropy increase, the solution above the diagonal line in Fig.|8]3 is eliminated when e up > 0.25 
while the solution below the diagonal line is eliminated when e up < 0.25. He also noticed that the jump behavior of 
an A-RD for e up > 0.25 is slow-mode-like (i.e. Sn and SB are anti-correlated) , while it is fast-mode-like (i.e. Sn and 
SB are correlated) for e up < 0.25. This directly relates to the fact that the jump of an A-RD equals the jump of the 
compound SS/RD wave (see Appendix IV. C), and a slow mode turns fast- mode- like when a > [4]. 



X. APPENDICES 



I - From anisotropic MHD to the Anisotropic DNLSB equation: 

From moment integrations of the Vlasov equation that neglect the off-diagonal components of the pressure tensor 
(the empirical validity of this approximation for our system is shown in Fig. 3(c) of Paper I), we can write down the 
anisotropic MHD (AMHD) equations p~5j[16]. The energy closure is undetermined. 

In Lagrangian form, 



d 
dt 



p + pd x V x = 



(26) 



P^V x +d x P + d x 

at 



Mo 



3 Mo 



(27) 



pi V t~— 9x(eB t ) = 
at po 



(28) 



d ( B 

— B t - B x d x V t + B t d x V x = d x (rj r d x B t ) - d x I — —e x x d x B t 

dt \Pone 



(29) 



d „ 7P d 



1 B 2 B 2 

-(s + 2)—-s^ 
3 /io Mo 



Mo 



d x V x + (7 - 1)(1 - e) 1 ^ ■ d x Y t + (7 - l)d x Q x = (30) 
Mo 
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where 



■2P± 



B 2 /Mo 



(31) 



7 = 5/3 or 7/5 for monatomic or diatomic plasma, respectively [32]. e, Pn, Pj_, p, T4, V^, 5^, B t , Q x and ?? r are 
the temperature anisotropy factor, pressure parallel to the local magnetic field, pressure perpendicular to the local 
magnetic field, mass density, velocity of the bulk flow in the normal direction (e x ), velocity of the bulk flow in the 
tangential direction (y-z plane), normal component of the magnetic field, tangential components of the magnetic field, 
the heat flux in the ^-direction and the magnetic resistivity (assumed constant). The first and the second term on 
the RHS of Eq. p9| ) are from the magnetic dissipation and the Hall term respectively. 

Then we follow the procedure of Kennel et al [4 , using Lagrangian mass spatial coordinates, 

-J d u x r =[ ^dx (32) 
at J p 

Jumping to the upstream (subscripted by "0") intermediate frame in order to separate the slow and fast variations, 
we take 

r y/e^C An t, V = x r - yfe^C An t (33) 
where C\ n = B 2 /(popo), and with the approximations 

dr^d^ 8b 2 t <^b 2 m ] fe<C£ ; ^P<Po; SP <C P , (34) 

(where 5 means variation), we can collapse the original seven equations into two coupled equations 

d T b t + d v [ab t (b 2 - b 2 t0 ) + nb t (e - e ) + Ab t (Q x - Q x0 )] = d^Rd^) - -^—d z d 2 (e x x b t ) (35) 



with 

,= N .. ^ „ 

12[(e - A)C\ n - C|] ' 6[(e -A)C 2 An - CI] ^ 2e ' 2^[(s - A)C\ n - C|] 



a _ [(3 7 -l) g0 -(37-4)]Ci„ ; (& f 2 o-2)Cj n | 1 ; A _ (7-l)Cj n (3g) 



and 

A=?( 7 -l)(l- e o)(l + J&) (37) 

where C% 7^b/Po, R Vr /(2y^C Un), and v / m i /(/i ^e 2 ). 

Since the heat flux is approximately proportional to V|P| (as pointed out in Fig. 3(a) of Paper I), it should 
enter the source term on the RHS as —Ad 2 \b t \. In plasmas with f) n > 1, A is negative and hence the heat flux helps 
shocks dissipate energy. We implicitly incorporate it into the resistivity R (as, for instance, is done for the shear and 
longitudinal viscosities discussed in [4 ). We then arrive at the anisotropic DNLSB equation, Eq. ([!]). This equation 
can also be derived from regular reductive perturbation methods with a proper ordering scheme. For instance, a 
DNLS equation with the CGL condition and more corrections, including finite ion Larmor radius effects and electron 
pressure, was derived using regular reductive perturbation methods [33] . 

II - Non-steepening and non-spreading of the intermediate mode: 

Beginning with the anisotropic MHD equations in Lagrangian form in Appendix I, we neglect the dissipation and 
the Hall term on the RHS of Eq. ( [29] ) . Then the simple wave solution can be obtained by substituting d/dt —> —C5, 
dx —> &x5: where C is the wave speed and S means variation. [34] 

- C5p + p5V x = (38) 
B 2 B 

- CpSV x + SP ± ^Se + —SB Z = (39) 

Mo Mo 

- CpSV z - e—5B z - -^Se = 40 
Mo Mo 
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Cp5Vy-e—8B y = 
Mo 



CSB Z + B Z SV X - B X 5V Z = 



(41) 
(42) 



C8B„ 



BJV y 







(43) 



Eq. (41) and Eq. (43) give us the intermediate speed Cj = \/sB x / JpoP- Combined with Eq. (38), the steepening 



tendency of an intermediate mode can then be expressed as, 



S(Ci + V x ) 



Cj_(5p Se 
2 \P e 



(44) 



Using Eqs. (38), (39), (40) and (42), we get 5{Ci + V x ) = 0. Therefore, the intermediate mode in anisotropic MHD 
does not steepen or spread, no matter what energy closure is used. It is linearly degenerate, as is its counterpart in 
isotropic MHD. 

Ill -The integral curves and Hugoniot Locus: 

To find the integral curves, we follow the eigenvector of the slow mode to form a curve, 



db z 



db y 



by, 



(45) 



where ( is a dummy variable. The integral curve is b y = (b y o/b z o)b z . 
For the intermediate mode, 



dby 
~d~C 



-b z . 



(46) 



-'to- 



Therefore, the integral curve is b 2 

As to the Hugoniot locus, we need to compute the shock speed S = (/i(q) — /i(qo))/(* — Qoi) where /, is the flux 
of Eq. ^ and i = 1 or 2. 



These can be combined to give, 



S(b z - b z0 ) = ab z (b 2 t - b 2 t0 ) + VLb z 5e 
S(b y - b y0 ) = ab y (b 2 - b 2 t0 ) + nbySe 

(b z0 b v - b y0 b z )[a(b 2 t - b 2 t0 ) + me] = 0. 



(47) 



(48) 



The first root is the Hugoniot locus of the slow mode: b y = (b y o/b z o)b z . For e(bt), so that Se ~ [d(Se)/d(b^)](b^ —b^ ), 
the second root gives us the Hugoniot locus of the intermediate mode: b\ = b^ . Although these results are the same 
as derived from the integral curves, this is not generally the case. 
IV - The pseudo-potential of Anisotropic MHD (AMHD): 

In the de Hoffmann- Teller frame, the jump conditions can be written as (following Hau and Sonnerup's procedure 



[pV x }° = 



(49) 



PV X 2 



£2 
Po 



3 

Mo 



(50) 



pV x V t 



B x ~B t 
Mo 







(51) 



^2 + ^_ p+ l (£ _ 1) ^, l , 

I 7-I 3 fio 



Mo Mo 



(52) 
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where we define a jump relation [Q]° = Qo — Q, with Q the upstream value and Q the value inside the transition 
region. From Eq. (49)- (52), we can derive 



A\ 



V 2 



* ~ BUM 



-b ± Vb 2 - 4ac 
2a 



(53) 



where 



a= 1 - 



7-1 

27 ' 



(54) 



2 



sec 6*0 + sq 



1 \ fi2 

2/ 



27-2 
37 



7 



1 , .B? 



(55) 



c = 



l-i S ec^o^o 
27 



go 2 7 -2 
2 37 



(£0 - 1) 



sec 2 #o^4xO 



7-1 
27 



(A 2 - £ )tan6*o + e 



B x 



(A^ -£ )tan(9 + £ 



(56) 



with cos#o = B x /B . 

The generalized Ohm's law is 



E + VxB=^J 



1 



-(JxB), 



(57) 



where the first and the second term on the RHS are the magnetic dissipation and the Hall term respectively. With 
the final jump condition [E t ]° = 0, we obtain 



4-^0(1 + h 2 )^ = (A 2 X - e)(B v - hB z ) + (A 2 x0 - s )(hB z0 - B y0 ) = F^amhd = 
tiQ ax oJd 



AMHD 



(58) 



^d i0 (l + h 2 )^ = (Al - e)(B z + hB v ) - (A 2 x0 - s )(hB y0 + B z0 ) = F Z , AUHD = ^amhd (59) 
rio ax Oid z 

where h = B x /ner] r measures the ratio of the dispersion to the resistivity. The pseudo-potential ^amhd is uniquely 
defined since dF ZiA MKD/dB y = dF yiA MHD/dB z . 

A: Fig. [9ja) shows the pseudo-potential of AMHD for the same parameters as Fig. [4ja) (which was calculated 
based on the reduced ADNLSB formulation). If Vs = then Al = so, and thus the potential minimum (where 
F Zt amhd = 0) occurs at A 2 = e. This implies that in the shock frame (also the upstream intermediate frame), 
A/, amhd = —V x + Cj = (— A x + y/e)B x / JJi^p = at the potential minimum. This is essentially the same point q m 
(where A/ = 0) of Fig.|4jd) with ADNLSB. Fig.[9jb) shows a similar reversal with the CGL closure. We note that the 
CGL closure exhibits an even stronger tendency to reverse the pseudo-potential. In Fig.[9^c), the pseudo-potential for 
Obn = 75°, /3o — 0.4 and C2 = 0.2 is shown, these parameters are more similar to those seen in our PIC simulation. 

B: Now we consider the jump conditions to an asymptotic downstream by neglecting the LHS and terms with h 



of Eqs. (58) and (|59|). The relation B t ^/B x — tan#o(A^ — sq)/(A 2 x d — Sd) can be derived where we label quantities 



Q ~~ ^ Qd ("d" for downstream). We can eventually invert A x0 as a function of A 2 d from Eq. (53). The result is 
plotted in Fig. [9td) which shows possible shock solutions as functions of the downstream intermediate Mach number 
Mj d = V 2 d /Cj d = A 2 d jsd [33]. In the green curve {e& = €o case), the portion from A-RD (anisotropic- RD) to SSS 
is the IS branch, from SSS to LS (linear slow mode) is the SS branch. When s& < £o> a new slow-shock transition 
from A x0 = So > 0.25 to A 2 X d = Sd is noted at the point (1,1). This new SS constitutes the slow shock portion of a 
compound SS/RD wave. For a given £q, the smallest possible £d, m in shrinks the SS and IS branches to the point (1,1). 
It can be shown that £d,min > 0.25 is always true for £q > 0.25. Therefore the existence of the new slow shock requires 
Sd > 0.25. In other words, £o > e m (= £d) > 0.25 is always true for a SS/RD compound wave in full anisotropic MHD. 

C: From Fig.[9^d) and further investigations, it can be shown that the anisotropic-RD(A-RD) at (1,1) has the same 
jump as that of the new SS at (1,1) plus a RD that does not change e and thermal states. Therefore an A-RD and 
the corresponding compound SS/RD wave have the same jump relations. 

V - The Equal- Area Rule and Intermediate Shocks: 
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The equal-area rule applies to conserved quantities in hyperbolic equations, which in our case is b z . From Eq. ([6| 
and the general form of Eq. (20), we find a simple relation between the pseudo- force and the flux function, 

F z \ by =o = -ah(b 2 z - b 2 z0 ) - nb z de(b z , b y ) + V s (b z - b z0 ) = -fi\ by=0 + V s (b z - b z0 ) (60) 
From Eq. (|6| and Eq. (J8|, a simple relation between the slow characteristic and the flux function is 

XsL ^= o = wX v =o> (61) 

It is then easy to show that, 

(Xs L \ by =o - V s )db z = -F z \ by=0 (62) 

This indicates that a stationary point b z , where F z = 0, will be located where the integral on the LHS is zero. This 
is called the equal-area rule. From this relation, with a given b z o and b Zl we can determine the shock speed Vs that 
causes the integral to vanish. Or for a given b z o and Vs, we can determine the possible downstream state b z . We apply 
it to the following examples to demonstrate the formation of intermediate shocks (which have a super-intermediate 
to sub-intermediate transition) in isotropic MHD. 

When the upstream q r = (b z o,0) is given and fixed, we can vary = (6^,0) to see the effect on possible shock 
solutions. In Fig. [lC^a), when is chosen above b z = 0, a slow shock solution is found by determining a proper 
horizontal line (q r — q/; note that the vertical position measures the shock speed Vs), which makes the red area below 
the line q r — q^ equal the red area above. The shock speed is slower than the upstream intermediate speed (black 
horizontal line across 0), the upstream (point q r ) is super-slow and sub-intermediate (Vs > Asx| qr , Vs < A/| qr . Since 
A equivalents to C — u where C is the phase speed and u is the bulk flow speed measured in upstream intermediate 
frame, Vs ^ A implies that the mach number measured in shock frame M = (Vs + u)/C ^ 1) and the downstream 
(point qj) is sub-slow (Vs < Asx| qz ). Traditionally in isotropic MHD, the super-fast state is termed number 1, sub-fast 
and super-intermediate is 2, sub-intermediate and super-slow is 3, and sub-slow is 4. Therefore a slow shock is also 
called a 3-4 SS. 



In Fig. 10 (b), if q/ is chosen below the point b z = 0, a 2-4 intermediate shock (q r — q^) is formed, with upstream 
being super-intermediate and downstream being sub-intermediate and sub-slow. In Fig. fToFc), with the same shock 
speed, a 2-3 IS transitions to a q/ with a more negative value is also possible. Note that the jump cross a compound 
2-3 IS/ 3-4 SS (from this q^ to the q/ in (b) ) equals to that of the 2-4 IS in (b). In Fig. [TO^d) , with the same q/ of 
Fig. [Tofc), a 2-3=4 IS (q r — q m ) with the maximum IS speed could be formed and attached by a slow rarefaction 
(q m — q/). This is a compound IS/SR wave, with 6 z | qrn = —b z o/2 which can also be determined by Asx(6 2 | qrr J = 
[fi(b z o) — /i(&s|q m )]/(fr*o — &z|q m )> as snown by Brio and Wu [2]. Similar arguments can be made in a system with 
fast and intermediate modes. 

Therefore, an intermediate shock is not directly associated with an intermediate mode. It is steepened by magneto- 
sonic waves (slow or fast modes), not by intermediate mode itself. This was first justified by Wu's (1987) coplanar 
simulations [3 (i.e., no out-of-plane magnetic field is allowed), where the intermediate shock forms even though the 
intermediate mode is not included (since the out-of-plane SB y is necessary for nontrivial solutions of the intermediate 
mode, as shown in Appendix II). The coupling of intermediate and magneto-sonic waves and the admissibility of 
intermediate shocks in the ideal MHD system was discussed by Kennel et al. [4]. 
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FIG. 1: [From Liu et a/., (2011) [J] The e distributions of runs with B n = 60° at 500/ft C2 , 75° at 200/0^, 83° at 700/0^. 
@bn is the angle between the magnetic field and the shock normal direction (e^) far upstream in these shock simulations. Q C i 
is the ion cyclotron frequency based on the upstream magnetic field. These simulations are designed to study the structure 
of reconnection exhausts in the normal direction. The downstream e tends to plateau at 0.25. When e < 0, the plasma is 
susceptible to the firehose instability. 




FIG. 2: Panel (a): An initial discontinuity between q r and results in two waves, the "1-wave" and "2- wave", that propagate 
in the r\ direction along time r. The middle state q m needs to be determined; Panel (b): The state space plot in the (b z ,b y ) 
plane. The value q r = (b z o, 0) is chosen since there is no out-of-plane B y upstream of the slow shocks in Paper I. q r straightly 
connects to qz and forms a switch-off slow shock; Panel (c): In order to connect q r to qz, it is necessary to cross the degeneracy 
band into the reversal region, which could cause the path to rotate at q m . 
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FIG. 3: The e distribution vs. 1/B 2 for the cases Obn = 30°(yellow), 45°(magenta), 52°(green), 60°(blue), 75°(red) and 
83° (black) from Fig. 5 of Paper I pQ. The dashed line has slope -0.5. In comparison, the diamond curve is the theoretical 
prediction with the CGL condition for the Obn — 75° case. 
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FIG. 4: Panel (a): A pseudo-potential ^ with Vs = 0. Upstream (point q r ), 0q = 42°, /3q = l,e = 1,C2 = 0.5 and we choose 
^(Qr) = 0. Since the transition occurs within the radius bt = b z o, we set ^ = for bt > b z o for a better visualization. The 
potential for negative B y is mirror symmetric to the part shown here; Panel (b): e(6t); Panel (c): Cuts of \P, F z , and a e ff along 
the B z axis with B y = 0; Panel (d): Aj and Xsl along the B z axis with B y = 0. The vertical axis measures speed (normalized 
to Cau)- "D" stands for degeneracy. The red area above Vs (zero here) equals the red area below Vs, and the same rule applies 
to the blue area. 




FIG. 6: The pseudo-potential with so = 0.24 < 0.25, while other parameters are the same as Fig.^] 
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£ B,s Hodogram 
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x/d ; x/d ; B z 

FIG. 7: Results of PIC simulations (runs f, g and h of Paper I 1 ]). Row (a): A case with Obn = 75°, fto = 0.4 and initial 
width Wi = ldi at time 200 /Q C i. e is shown on the left, different magnetic components in the middle, B z -B y hodogram on the 
right; Row (b): A similar case with a wider initial width Wi = 10di at time 450/^ C i; Row (c): A similar case to (a), but with a 
weak guide field B y o = 0.2Bo at time 200 /Q C i- The dotted curves in the center column are the predicted B z magnitudes and 
positions of switch-off slow shocks (SSS) or rotational discontinuities (RD) from isotropic MHD theory [31] . 
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FIG. 8: In panel (a), the possible jumps of e of an anisotropic- RD, are constrained by requiring positive P±d, P\\d an d B 2 d . The 
red region is forbidden. In panel (b), the plot is further constrained by the requirement that entropy increases (with an entropy 
from the H-theorem defined as ln(P|^ 3 P^ 3 /p 5 / 3 ) for a bi-Maxwellian distribution). The constraint of increasing entropy has 
eliminated the region above the diagonal line when s up > 0.25, and the region below the diagonal line when e up < 0.25. Here 
/3q = 1; a higher /3q would collapse the valid region into a narrower region along the diagonal line. 
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FIG. 9: Results with fully anisotropic MHD. Panel (a) is the pseudo-potential ^ with our closure, Eq. (15). Upstream (point 
qr), Vs — 0, 0o = 42°, /3q = 1, so = 1 and C2 = 0.5; Panel (b) is ^ with the CGL closure. Upstreani(point q r ), Vs — 0, 
6> = 42°, f3 = 1 and e = 1.5; Panel (c) is ^ for = 0, 6 = 75°, /3 = 0.4, e = 1, and c 2 = 0.2 with our closure; Panel 
(d) is the shock curve with upstream parameters 0q — 42°, f3 = 1 and So = 1. In the green curve (s^ = So case), the portion 
from A-RD (anisotropic-RD) to SSS is the IS branch, from SSS to LS (linear slow mode) is the SS branch. Different curves 
represent cases with different Ed of values 1, 0.95, 0.9, 0.85 (from outer curve to inner curve). Other than the A-RD, a new SS 
exists at (1,1) when Sd < £o- Both the IS and SS branches shrink toward the point (1,1) as Sd decreases. 




FIG. 10: Application of equal-area rules with cases in ideal (isotropic) MHD. The A/ and Xsl along b z are measured in the 
upstream intermediate frame with upstream (point q r ) parameters, 0o = 42°, f3o = l,£o = 1. The vertical axis measures speed 
(normalized to Cau)- Once the qz is chosen, the shock speed Vs (measured by the red horizontal line) can be determined by 
equating area (between the Xsl and the red line) above the red line to area below the red line. This rule results in shock speeds 
(a) V s = -0.1; (b) V s = 0.05; (c) V s = 0.05; (d) V s = 0.0958. 



